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A  planar  slab  of  negative-index  material  works  as  a  superlens  with  sub-diffraction-limited 


resolution,  as  propagating  waves  are  focused  and,  moreover,  evanescent  waves  are 
reconstructed  in  the  image  plane.  Here  we  demonstrate  a  superlens  for  electric  evanescent 
fields  with  low  losses  using  perovskites  in  the  mid-infrared  regime.  The  combination  of  near¬ 
field  microscopy  with  a  tunable  free-electron  laser  allows  us  to  address  precisely  the  polariton 
modes,  which  are  critical  for  super-resolution  imaging.  We  spectrally  study  the  lateral  and 
vertical  distributions  of  evanescent  waves  around  the  image  plane  of  such  a  lens,  and  achieve 
imaging  resolution  of  A/14  at  the  superlensing  wavelength.  Interestingly,  at  certain  distances 
between  the  probe  and  sample  surface,  we  observe  a  maximum  of  these  evanescent  fields. 
Comparisons  with  numerical  simulations  indicate  that  this  maximum  originates  from  an 
enhanced  coupling  between  probe  and  object,  which  might  be  applicable  for  multifunctional 
circuits,  infrared  spectroscopy  and  thermal  sensors. 
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In  1968,  Veselago  theoretically  discussed  a  material  with  a  negative 
refractive  index,  as  well  as  its  extraordinary  responses  to  electro¬ 
magnetic  waves,  including  negative  refraction,  reversed  Doppler  and 
Cerenkov  effects  and  focusing  with  a  planar  lens1.  Pendry  revisited  this 
idea  30  years  later,  showing  that  with  such  a  planar  lens  one  can  create 
an  image  with  a  resolution  beyond  the  conventional  diffraction  limit2. 
This  opened  the  rapidly  growing  field  of  metamaterials,  which  present 
opportunities  for  new  and  remarkable  applications  in  electromagnet¬ 
ics3,4  ranging  from  electrostatics5,6  via  radio  frequency7,  microwave8 
and  terahertz9  to  the  infrared10  and  optical  regime11-14.  The  experi¬ 
mental  proofs  of  negative  refraction  and  sub-diffraction-limited  reso¬ 
lution  by  a  negative-index  planar  lens  were  first  reported  in  the  GHz 
range15,16.  Subsequently,  it  was  shown  for  visible17  and  infrared18  wave¬ 
lengths  that  the  evanescent  field  information  of  transverse-magnetic 
(TM)  waves  can  be  recovered  using  a  medium  with  negative  permitti¬ 
vity  (e)  only.  Such  a  lens  is  normally  called  a  superlens,  and  substantial 
progress  has  been  achieved  in  this  rapidly  developing  area19. 

A  superlens  is  preferably  realized  by  low-loss  materials2,19, 
among  which  perovskite  oxides  are  good  candidates.  In  the  mid- 
infrared  range,  perovskite  oxides  show  phonon  resonances,  which 
are  located  at  slightly  different  wavelengths  for  related  materials 
such  as,  for  example,  bismuth  ferrite  (BiFe03)  and  strontium  titan- 
ate  (SrTi03)20,21.  On  the  high-frequency  side  of  these  phonon  reso¬ 
nances,  perovskites  exhibit  negative  permittivities,  which  are  suit¬ 
able  for  superlens  structures2,19.  At  these  wavelengths,  the  intrinsic 
absorption  of  light  is  small20,21.  In  addition,  epitaxially  grown  oxides 
exhibit  highly  crystalline  interfaces  resulting  in  low  scattering.  Both 
the  small  intrinsic  and  low  scattering  losses  of  perovskites  could  sig¬ 
nificantly  improve  the  imaging  resolution  of  superlenses. 

Perovskite  oxides  exhibit  many  intriguing  properties  such  as 
colossal  magnetoresistance22,  ferroelectricity23,  superconductivity24 
and  spin-dependent  transport25,  which  lead  to  numerous  applica¬ 
tions  such  as  non-volatile  memories26,  microsensors  and  microactu¬ 
ators27,  as  well  as  in  nanoelectronics28.  Piezoelectricity  and  ferroelec¬ 
tricity  allow  for  the  manipulation  of  geometrical  dimension,  electric 
polarization  and  dielectric  properties  by  external  electric  fields29. 
Multiferroic  perovskites,  such  as  BiFe03,  exhibit  ferroelectricity  as 
well  as  magnetoelectric  coupling29,  thus  providing  a  pathway  for 
additional  degrees  of  tunability.  Moreover,  perovskites  with  match¬ 
ing  lattice  constants  can  be  grown  epitaxially  on  top  of  each  other30,31, 
which  allows  one  to  combine  their  properties  in  multifunctional  het- 
erostructures31-33.  These  unique  properties  of  perovskites  may  inspire 
new  functionalities  of  metamaterial  structures  such  as  tunability  by 
external  fields  and  novel  designs  for  multifunctional  circuits. 

Here,  we  study  a  new  type  of  superlens  for  electric  fields  in  the 
mid-infrared  based  on  the  perovskite  oxides  BiFe03  and  SrTi03.  We 
investigate  the  evanescent  waves  in  the  image  plane  of  perovskite 
superlenses  by  means  of  scattering-type  near-field  infrared  micro¬ 
scopy  (s-NSIM)34,35.  The  combination  with  a  free-electron  laser 
(FEL),  which  is  precisely  tunable  in  the  wavelength  (A)  regime  from 
4  to  250  jam,  enables  us  to  address  the  polariton  modes  and  study 
their  lateral,  vertical  and  spectral  distribution.  Such  detailed  char¬ 
acterizations  are  important  to  understand  superlensing  effect,  but 
they  were  not  comprehensively  conducted  in  previous  work17,18.  At 
certain  wavelengths  we  observe  enhanced  evanescent  fields  in  the 
image  plane  as  well  as  a  resolution  beyond  the  classical  diffraction 
limit.  Moreover,  we  find  that  the  evanescent  fields  show  maxima  at  a 
certain  distance  between  the  probe  and  the  sample  surface,  for  which 
the  vertical  position  depends  on  the  wavelength.  Comparisons  with 
two-dimensional  numerical  simulations  indicate  that  a  superlens- 
enhanced  coupling  between  probe  and  object  causes  this  effect, 
which  might  be  utilizable  for  controlling  polariton  propagation. 

Results 

Perovskite-based  superlens.  The  original  superlens  proposed  by 
Pendry  is  a  single  slab  of  a  material  with  £=  - 1  surrounded  by 


Figure  1 1  s-NSIM  setup  and  perovskite  properties,  (a)  Sketch  of  the 
experimental  setup  including  the  superlens,  the  geometry  at  the  near¬ 
field  probe  (blue),  and  the  free-electron  laser  light  source.  The  superlens 
consists  of  the  layers  A  (BiFe03)  and  B  (SrTiOs)  of  thicknesses  d  and  2d 
(d  =  200  nm),  respectively.  The  objects  to  be  imaged  are  SrRu03  patterns 
(orange  rectangles)  on  a  SrTi03  substrate.  All  constituents  of  the  superlens 
are  perovskite  oxides  that  match  in  their  crystalline  structures  resulting 
in  low  scattering  at  the  highly  crystalline  interfaces.  The  near-field  tip 
probes  the  evanescent  fields  on  the  image  side  of  the  lens.  The  superlens 
is  excited  by  an  infrared  free-electron  laser,  which  is  precisely  tunable  in 
the  range  from  4  to  250  pm.  (b)  Imaginary  and  real  parts  of  the  dielectric 
constants  eof  all  constituents  (SrRuOs  data  determined  by  Fourier 
transform  infrared  spectroscopy,  values  for  SrTiOs  and  BiFeOs  taken  from 
literature2021,  see  Supplementary  Fig.  SI);  (c)  real  parts  of  the  dielectric 
constants  at  the  high-frequency  side  of  their  phonon  resonances  depicted 
in  b.  The  arrow  indicates  the  wavelength  at  which  superlensing  is  expected. 


air  (£air=  +  1)2.  At  the  two  interfaces  of  the  slab,  coupled  surface 
polariton  modes  are  excited,  which  amplify  the  evanescent  electric 
fields  arising  from  an  object  and  transform  them  to  the  opposite 
side  of  the  lens2.  The  oxide-based  superlenses  in  our  study  consist  of 
layers  of  matched  perovskites  (see  Fig.  la).  We  structure  a  50-nm- 
thick  film  of  metallic  perovskite  strontium  ruthenate  (SrRu03)  on  a 
SrTi03  substrate  by  photo -lithography,  acting  as  superlens  objects 
with  a  lateral  size  down  to  3  pm.  To  ensure  a  fixed  distance  between 
obj  ects  and  the  slab  we  add  a  layer  A  (thickness  d ,  £A=£air  =  +  1 )  on  top 
of  the  objects  before  we  grow  the  actual  superlens  layer  with  £B=-  1 
(thickness  2d).  As  layer  A  has  a  slightly  different  permittivity  than 
air  this  two-layer  superlens  is  asymmetric.  Hence,  we  also  study 
a  superlens  with  a  symmetric  design  after  adding  an  additional 
layer  of  material  A  (thickness  d)  on  top.  Superlensing  is  expected 
for  such  systems,  when  the  real  parts  of  the  permittivities  DAc(£a) 
and  DAe(£B)  have  the  same  absolute  values  and  opposite  signs 
93c(£b)  =  -  9\e(eA)  <  0.  Assuming  the  dielectric  constants  of  the  layers 
are  similar  to  the  ones  described  in  literature20,21  (see  Supplementary 
Fig.  SI),  this  superlensing  condition  is  fulfilled  at  A=  13.9pm  with 
the  layers  A  and  B  being  BiFe03  and  SrTi03,  respectively  (Fig.  lb,c). 
Many  other  similarly  suitable  pairs  of  perovskites,  such  as  PbZrJq, 
x03  and  SrTi03,  BiFe03  and  PbZrxTi1_x03,  BaTi03  and  SrTi03,  or 
BiFe03  and  BaTi03,  can  be  found  at  different  wavelengths,  implying 
the  flexibility  of  perovskite  superlenses  in  terms  of  the  operation 
wavelength.  Please  note  that  other  groups  of  dielectrics  such  as,  for 
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Figure  2  |  Near-field  images  of  three  different  samples.  From  left  to  right,  the  subfigures  display:  sketches  of  the  sample,  topography  images  obtained  by 
atomic-force  microscopy  (scale  bars,  10pm),  near-field  signals  as  functions  of  the  probe-sample  distance  z  for  selected  wavelengths  (an  offset  is  added 
for  better  comparison),  as  well  as  near-field  images  as  described  in  the  following,  (a)  For  SrRuOs  objects  (orange  rectangle)  on  a  SrTi03  substrate  we 
image  second-  and  third-harmonic  near-field  signals  (NF2Q  and  NF3Q)  using  a  C02  laser  (A  =  10.6pm).  (b,c)  For  both  types  of  superlenses  we  depict  NF3Q  at 
two  different  wavelengths  being  A= 17.5  and  14.6 Jim  for  the  symmetric  superlens  (b)  and  A= 17.3  and  14.1  jam  for  the  asymmetric  lens  (c).  The  red  and  green 
curves  in  the  distance  curves  correspond  to  areas  with  and  without  SrRu03  objects  on  the  opposite  side  of  the  lens,  respectively.  Such  distance  curves  and 
the  near-field  images  show  strong  signals  at  both  wavelengths,  but  only  at  the  shorter  wavelengths  we  observe  a  contrast  beyond  the  diffraction  limit  due 
to  the  superlensing  effect. 


example,  fluorides36,37  and  simple  oxides38  can  form  similar  pairs  for 
superlenses,  which  altogether  cover  an  even  larger  wavelength  range 
of  operation. 

We  examine  the  enhanced  evanescent  fields  on  the  image  side 
of  the  lens  with  s-NSIM  in  combination  with  a  FEL  (see  Fig.  la  and 
Methods),  which  allows  for  polariton-enhanced  s-NSIM39,40.  Sur¬ 
face  polariton  modes  can  be  excited  at  the  interface  of  two  materi¬ 
als,  which  have  permittivities  of  opposite  signs41,42.  In  the  studied 
wavelength  range  from  13.5  to  18.5  pm,  the  layers  A  and  B  of  the 
perovskite-based  superlenses  show  the  following  permittivities  (see 
Fig.  lb,c)  and  corresponding  polariton  modes:  when  A>  16.5  pm  we 
find  DT(£B)  <  9T(£A)  <  0  and  polariton  modes  are  excited  at  the  entire 
sample  surface,  being  the  interface  between  the  toplayer  and  air.  For 
A  <  16.5  pm,  it  is  9T(£B)  <  0  <  9A(£a)  and  polariton  modes  are  created 
at  both  interfaces  of  layer  B.  As  layer  B  is  only  400-nm  thick,  these 
two  modes  are  coupled,  but  are  not  necessarily  localized.  We  observe 
an  enhanced  near-field  signal,  whenever  a  polariton  mode  is  present 
at  the  position  of  the  probe.  However,  only  around  the  superlens¬ 
ing  wavelength  (£A=  -  £B)  the  polartion  modes  are  strongly  localized 
and,  hence,  create  a  sub -diffraction -limited  image  of  the  objects  on 
the  opposite  side  of  the  lens.  Please  note  that  with  near-field  image, 
we  refer  to  the  image  obtained  by  near-field  microscopy.  The  con¬ 
trast  in  these  images  arises  from  local  sample  properties  and  thus 
different  coupling  to  the  probe,  allowing  us  to  distinguish  the 
different  objects  on  the  opposite  side  of  the  superlens  with  a  sub- 
diffraction -limited  resolution  (see  Methods)43. 

Mid-infrared  near-field  imaging.  When  placing  the  probe  at  a  dis¬ 
tance  of  about  30  nm  to  the  sample,  which  is  scanned  relatively  to 


the  probe  position,  we  obtain  the  near-field  response  as  well  as  the 
topography  of  the  sample.  The  corresponding  results  are  depicted 
in  Figure  2  for  similar  object  distributions  with  lateral  sizes  from 
3x3  to  8x8  pm2  and  three  different  samples,  namely  (a)  structured 
SrRu03  objects  on  a  SrTi03  substrate,  (b)  a  symmetric  superlens 
and  (c)  an  asymmetric  superlens  without  the  top  BiFe03  layer. 

The  SrRu03  objects  without  a  lens  (Fig.  2a)  show  a  near-field 
response,  which  depends  on  the  polarization  and  the  incident  angle 
of  the  impinging  light  with  respect  to  the  objects.  The  k-vector  direc¬ 
tion  of  the  incident  TM-polarized  light  is  shown  in  the  sketch  and  in 
the  topography  image  in  (a)  and  is  the  same  for  all  scans  depicted  in 
Figure  2.  Over  a  wide  wavelength  range  in  the  infrared,  the  permit¬ 
tivity  of  SrRu03  is  much  smaller  than  zero,  resulting  in  an  excellent, 
metal-like  scattering  behaviour.  Consequently,  the  distance  depend¬ 
ence  of  the  near-field  signal  on  SrRu03  exhibits  the  same  exponen¬ 
tially  decaying  character  of  the  evanescent  waves  for  all  wavelengths 
(see  middle  panel  of  Fig.  2a).  The  metallic  SrRu03  structures  show 
enhanced  field  intensities  when  their  geometries  match  the  wave¬ 
length,  an  effect  which  is  well  known  in  plasmonics44.  Consequently, 
certain  parts  of  the  objects  appear  bright  in  the  second-  and  third- 
harmonic  near-field  signals  NF2Q  and  NF3Q  (see  Methods).  In  gen¬ 
eral  NF2Q  is  about  three  times  larger  than  NF30>  with  an  interference - 
like  background  signal.  To  ensure  pure  near-field  detection  in  the 
following,  we  will  depict  third-harmonic  signals  only. 

For  the  symmetric  superlens  (Fig.  2b),  at  A=  17.5  pm  we  observe 
a  near-field  signal  due  to  a  non-localized  polariton  mode  at  the  sam¬ 
ple  surface,  showing  no  clear  contrast  between  areas  with  and  with¬ 
out  SrRu03  objects.  Close  to  the  superlensing  wavelength,  at  around 
14.6  pm,  the  evanescent  fields  are  highly  localized,  resulting  in  a 
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Figure  3  |  Spectral  response  of  the  symmetric  superlens.  Near-field  spectra  as  well  as  near-field  images  and  NF3Q-distance  curves  for  selected  wavelengths 
(scale  bars,  10pm).  (a)  For  A=  13.9-15.9  Jim,  an  imaging  contrast  exists  because  of  the  localized  polariton  mode,  (b)  Shows  the  near-field  spectrum  for  a  fixed 
distance  of  z  =  20nm  with  the  results  between  13.5  and  16.25  pm  being  multiplied  by  a  factor  of  four  in  the  plot,  (c)  For  A= 16.8-18.4  jim,  no  imaging  contrast 
is  observed  although  the  near-field  signal  is  enhanced  by  the  non-localized  polariton  mode.  The  red  and  green  curves  in  the  distance  curves  of  a,  c  and  the 
spectrum  b  correspond  to  the  near-field  signals  on  areas  with  and  without  SrRuOs  objects  on  the  opposite  side  of  the  lens,  respectively. 


strong  contrast  reconstructing  the  SrRu03  structures.  Please  note  that 
even  though  the  NSIM  probe  scans  more  than  800  nm  higher  above 
the  SrRu03  objects  compared  with  the  bare  SrRu03  case,  the  signal  is 
about  two  times  stronger  because  the  coupled  polariton  modes  at  the 
interfaces  enhance  the  near-field  signals  arising  from  the  objects. 

On  the  asymmetric  superlens  (Fig.  2c),  we  observe  a  similar, 
slightly  blue-shifted  response  compared  with  the  symmetric  super¬ 
lens,  with  non-localized  polariton  modes  around  17.3  jam  and 
superlensing  with  maximum  contrast  at  around  14.1  pm.  The  signal 
on  the  asymmetric  superlens  is  two  times  stronger  than  for  the  sym¬ 
metric  superlens,  as  the  near-field  probe  is  placed  by  200  nm  closer 
to  the  top  SrTi03  surface  at  which  the  polariton  mode  is  excited. 

With  both  superlenses  we  clearly  resolve  the  buried  SrRu03 
objects  at  the  corresponding  superlensing  wavelengths.  In  the  fol¬ 
lowing,  we  discuss  in  detail  the  spectral  response  of  both  samples 
over  a  broad  wavelength  range. 

Spectral  response  of  the  near-field  signals.  Our  s-SNIM  setup 
allows  for  studying  polariton  modes  with  different  characteristics  in 
the  wavelength  range  from  13.5  to  18.4  pm.  For  the  symmetric  super¬ 
lens  (Fig.  3),  we  observe  an  enhanced  localized  signal  for  A  from  13.9 
to  15.5  pm  (Fig.  3a),  close  to  the  superlensing  condition.  The  high¬ 
est  contrast  is  observed  at  14.6  pm  and  the  signal  drops  to  zero  for 
wavelength  1  pm  smaller  or  larger  than  that.  Note  that  areas  without 
objects  show  a  small  negative  signal  in  the  superlensing  regime,  pos¬ 
sibly  due  to  evanescent  waves  on  the  substrate- superlens  interface 
with  opposite  phase.  At  longer  wavelengths,  both  layers,  BiFe03  and 
SrTi03,  show  negative  permittivities  and  a  single  polariton  mode  is 
excited  at  the  top-most  sample  surface  leading  to  an  enhanced  near- 
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field  signal  on  the  entire  sample  area  for  A>  17  pm.  Even  though  this 
signal  is  up  to  eight  times  higher  than  the  superlensed  signal,  it  does 
not  show  any  image  contrast  except  for  some  artifacts  at  the  topogra¬ 
phy  edges  caused  by  scattering  and  errors  in  the  distance  control. 

For  the  asymmetric  superlenses  (Fig.  4),  we  observe  a  local¬ 
ized  near-field  signal  around  the  superlensing  wavelength  for 
13.7  pm  <  A<  15.2  pm  with  a  clear  object-related  contrast  at  14.1  pm. 
In  addition,  the  top-most  layer  SrTi03  of  this  lens  supports  a  non- 
localized  polariton  mode  for  14 pm <  A<  15.5 pm.  In  contrast  to  the 
symmetric  case,  in  which  this  mode  is  spectrally  well  separated  from 
the  superlensing  wavelength,  for  the  asymmetric  lens  the  wave¬ 
length  regimes  of  superlensing  and  non-localized  polariton  overlap. 
Interestingly  both  modes  respond  with  opposite  phases  resulting  in 
a  destructive  superposition  on  the  objects,  which  appear  dark  in  the 
range  from  14.5  to  15.2  pm.  With  increasing  wavelength  the  non- 
localized  mode  dominates  and  the  images  are  blurred. 

The  smallest  structures  resolved  on  the  superlenses  are  about 
3x3  pm2  in  size  and  about  1pm  apart  corresponding  to  a  relative 
resolution  of  about  A/ 14.  This  sub -diffraction-limited  image  is 
created  by  the  superlensing  effect.  Compared  with  the  image  of 
bare  SrRu03  objects  (Fig.  2a),  the  superlensed  image  appears  homo¬ 
geneously  bright  as  the  incident  light  is  scattered  by  the  probe  and 
from  there  directed  towards  the  object  (see  Methods).  In  general, 
the  imaging  with  s-NSIM  is  influenced  by  the  probe  and  its 
vertical  position  as  it  is  in  particular  observed  in  the  cross-sections 
of  Figure  5,  which  will  be  discussed  in  the  following  paragraphs. 

Mid- infrared  near- field  cross-sections  and  spectroscopy.  We  focus 
on  the  asymmetric  superlens  to  study  the  vertical  distributions  of 
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Figure  4  |  Spectral  response  of  the  asymmetric  superlens.  Near-field  spectra  as  well  as  near-field  images  and  NF3Q-distance  curves  for  selected 
wavelengths  (scale  bars, 10pm).  (a)  For  X-  13.7-15.2  jam,  an  imaging  contrast  exists  because  of  the  localized  polariton  mode,  (b)  Shows  the  near-field 
spectrum  for  a  fixed  distance  of  z  =  20  nm.  (c)  For  X-  15.4-18.4  jam,  no  imaging  contrast  is  observed  although  the  near-field  signal  is  enhanced  by  the 
non-localized  polariton  mode.  The  red  and  green  curves  in  the  distance  curves  of  a,  c  and  the  spectrum  b  correspond  to  the  near-field  signals  on  areas 
with  and  without  SrRuOs  objects  on  the  opposite  side  of  the  lens,  respectively. 


evanescent  waves.  Similar  measurements  for  the  symmetric  super¬ 
lens  can  be  found  in  the  Supplementary  Figure  S2.  Figure  5a  depicts 
normalized  vertical  cross-sections  on  the  image  side  of  the  sample 
for  different  wavelengths  showing  the  near-held  signal  as  a  func¬ 
tion  of  the  probe-sample  distance  z  and  the  position  of  the  sample. 
The  topography  of  the  sample  is  reflected  by  the  dark  area  on  the 
bottom  of  the  pictures  with  a  4-pm  wide  SrRu03  object  in  the  cen¬ 
tre.  For  X>  14.8  jam,  the  probe  excites  propagating  polariton  modes 
close  to  the  SrTi03  surface.  Around  the  superlensing  wavelength, 
for  X-  14.5-13.5  pm,  the  evanescent  held  is  localized  on  the  object 
with  decreasing  signal  and  contrast  for  shorter  wavelengths.  For  the 
latter  wavelength  regime,  we  observe  an  intriguing  phenomenon: 
a  maximum  in  the  evanescent  held  appears  at  a  certain  distance  z0 
between  tip  and  sample  surface.  Moreover,  z0  increases  with  smaller 
wavelength  and  appears  at  distances  of  up  to  150  nm  from  the  sam¬ 
ple  surface  at  X-  13.5  pan.  At  hrst  glance,  this  effect  seems  to  be 
unexpected  because  the  phonon-polariton  mode  is  a  conhned  sur¬ 
face  mode,  which  exponentially  decays  from  the  interface.  To  make 
sure  that  this  effect  is  not  an  artifact  in  our  experimental  setup,  we 
compare  these  results  with  numerical  simulations  as  discussed  in 
the  following. 

Figure  5b  shows  the  simulated  data  for  a  superlens  consisting 
of  the  same  constituents  and  geometry  as  in  our  experiments  (see 
Methods).  In  these  simulations  we  clearly  observe  the  same  effect  as 
in  Figure  5a,  that  is,  the  held  maximum  gradually  shifts  away  from 
the  sample  surface  when  the  wavelength  decreases.  This  maximum 
has  an  asymmetric  shape,  locating  on  the  right-hand  side  of  the 
structure,  which  is  illuminated  from  the  left  with  an  incident  angle 


of  75°.  The  asymmetry  might  be  formed  by  shadowing  of  the  struc¬ 
ture  by  the  probe:  when  the  probe  is  placed  on  the  left  side  of  the 
structure,  it  reflects  the  incident  light  and  light  hardly  reaches  the 
full  structure;  whereas  when  the  probe  is  positioned  on  the  right, 
the  structure  is  completely  illuminated  by  the  beam.  Compared 
with  the  experimental  results,  the  simulations  show  some  differ¬ 
ences:  First,  the  wavelengths  as  well  as  the  z-position  of  the  maxi¬ 
mum  are  slightly  shifted.  This  is  likely  due  to  the  small  discrepancy 
between  the  dielectric  constants  of  the  fabricated  superlens  layers 
and  the  data  reported  in  the  published  literature20,21,  which  are 
used  for  the  simulations  (see  Supplementary  Fig.  SI).  Second,  the 
near-field  signal  above  the  topographic  step  for  X>  14.5  pm  in  the 
simulation  shows  a  saddle  shape.  The  enhanced  signal  at  the  topo¬ 
graphic  steps  is  possibly  caused  by  the  sharp  edges  in  the  modelled 
geometry,  which  give  rise  to  a  highly  nonlinear  near-held  because 
of  localized  corner  polariton  modes45.  Please  note  that  this  effect 
is  only  relevant  in  the  propagating  polariton  regime,  whereas  it  is 
much  less  pronounced  around  the  superlensing  wavelengths  as  the 
polariton  mode  at  the  hat  interface  of  the  step  becomes  also  highly 
localized.  In  the  Supplementary  Figure  S3,  we  compare  these  results 
with  corresponding  simulations  without  the  topographic  step.  We 
hnd  very  similar  spectral,  lateral  and  vertical  response,  but  the 
enhanced  signal  at  the  edges  disappears.  This  comparison  verihes 
that  the  small  topography  on  the  sample  interfaces  because  of  the 
sample  fabrication  has  negligible  inhuence  on  the  observed  optical 
signals.  The  enhanced  coupling  of  probe  and  object  by  the  superlens 
layer  is  indeed  correlated  to  the  material  property  rather  than  the 
topological  protrusion. 
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Figure  5  |  Normalized  cross-sections  on  an  asymmetric  superlens,  (a)  Depicts  experimental  data,  and  (b)  shows  the  results  from  numerical  simulations 
(for  details  see  text).  The  horizontal  range  in  all  cross-sections  is  12pm  (scale  bars,  6jim).  The  topography  of  the  sample  and  the  position  of  the  SrRuOs 
object  is  reflected  by  the  dark  areas  at  the  bottom  of  the  figures.  The  localized  evanescent  fields  on  the  objects  show  a  maximum  at  a  certain  distance  z0, 
which  increases  with  smaller  wavelength.  This  effect  is  observed  experimentally  as  well  as  in  the  simulations  and  corresponds  to  a  superlens-enhanced 
coupling  of  probe  and  object.  (c,d)  Spectral  behaviour  of  the  near-field  maxima  in  the  experiment  (c)  and  in  the  simulations  (d)  for  areas  with  (red)  and 
without  (green)  objects  as  well  as  the  corresponding  contrast  V  (yellow)  calculated  from  this  data.  The  marked  areas  correspond  to  the  observation  of 
non-localized  near-field  signals  (green),  and  localized  evanescent  fields  due  to  the  superlensing  effect  (red). 


We  further  compare  the  results  of  experiment  (Fig.  5c)  and  simu¬ 
lation  (Fig.  5d)  by  analysing  the  absolute  values  of  the  maximum 
near-field  signals  with  and  without  SrRu03  object,  NFW  and  NFw/o, 
respectively,  as  well  as  the  resulting  contrast  V=  (NFw-NFw/o)/ 
(NFw  +  NFw/o).  Without  object  (green  curves)  we  observe  maxima  at 
15.3  pm  (c)  and  14.5  pm  (d)  due  to  non-localized  polariton  modes, 
whereas  on  the  SrRu03  objects  (red)  an  enhanced  signal  is  observed 
for  slightly  shorter  A.  The  resulting  contrasts  (yellow)  show  maxima 
at  14  pm  (c)  and  13.3  pm  (d),  respectively,  which  is  in  good  agree¬ 
ment  with  the  predicted  superlensing  wavelength  of  13.9  pm. 

Numerical  simulations  of  superlens-enhanced  coupling.  Why 

do  we  observe  an  enhanced  signal  at  a  certain  distance  to  the  sam¬ 
ple  surface?  The  coupled  phonon-polariton  modes,  which  create 
the  sub -diffraction -limited  image,  result  in  a  field  which  decreases 
exponentially  with  the  distance.  However,  as  we  place  a  scattering 
probe  on  the  image  side,  this  probe  itself  acts  as  an  object  next  to  the 
superlens,  creating  additional  fields  at  the  position  of  the  SrRu03 
objects.  The  probe-sample  system  therefore  consists  of  two  coupled 
scatterers  with  a  superlens  structure  between  both  of  them,  showing 
a  resonance  as  a  function  of  the  wavelength  and  the  probe-sample 
distance  as  discussed  in  the  following  paragraphs. 

Figure  6a  plots  the  electric-field  distribution  of  an  asymmetric 
BiFe03-SrTi03  superlens  with  a  line  source  in  the  object  plane  for 
different  wavelengths  and  no  probe  on  the  image  side.  For  wave¬ 
lengths  larger  and  smaller  than  the  superlensing  wavelength  around 
A  =  13.5  pm,  we  observe  unconfined  evanescent  waves  on  the  image 
side  of  the  superlens  due  to  the  excitation  of  non-localized  surface 
polariton  modes  at  the  SrTi03-air  interface.  For  A=  13.5  and  14  pm, 
the  superlensing  effect  takes  place  with  an  enhanced  confined  field 
on  the  image  side  of  the  lens  and  localized  polariton  modes  at 
both  interfaces  of  the  SrTi03  layer.  The  superlensing  effect  is  fur- 
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ther  verified  by  the  transfer  function  simulation  of  electric  fields. 
Figure  6b  shows  the  isothermal  contour  of  the  transfer  function 
of  the  asymmetric  superlens  (for  the  transfer  function  of  the  sym¬ 
metric  superlens  see  Supplementary  Fig.  S4).  Around  13.5  pm,  the 
amplitude  of  TM  waves  at  the  top  interface  of  the  superlens  still 
maintains  reasonably  large,  even  for  a  tangential  wave  vector  up  to 
10  k0  (Fig.  6c).  In  contrast,  the  field  intensity  dramatically  decays  for 
large  wave  vectors  in  a  control  sample,  in  which  the  400  nm  SrTi03 
film  is  replaced  by  a  400  nm  BiFe03  layer. 

Figure  6a  shows  that  the  electric  field  at  the  SrTi03-air  interface 
possesses  in  all  cases  a  maximum  at  the  sample  surface,  reflecting  the 
exponentially  decaying  behaviour  of  the  polariton  mode.  The  situa¬ 
tion  changes  with  the  presence  of  a  scattering  probe  on  the  imaging 
side  of  the  sample.  Figure  6d  shows  a  simulation  for  the  ideal  case  of 
two  equally  sized  metal  spheres,  object  and  probe,  separated  by  the 
two-layer  superlens  for  a  fixed  wavelength  of  14  pm  and  for  different 
gaps  z  between  the  top  scatterer  (that  is,  the  probe)  and  the  sample 
surface.  For  small  gaps  z=25nm,  one  observes  a  large  field  under¬ 
neath  the  probe-scatterer  £tip  due  to  near-field  enhancement  as  well 
as  a  localized  polariton  mode  Eint  at  the  interface  between  SrTi03 
and  BiFe03.  When  the  gap  is  increased  to  75  nm,  Eti  decreases, 
but  Eint  is  much  larger  than  that  for  z  =  25nm,  whereas  both  fields 
decrease  with  z  for  gaps  larger  than  100  nm. 

In  Figure  6e,  we  plot  the  z-dependence  of  Eti  (green)  and  Eint 
(red),  as  well  as  of  the  integrated  Poynting  vector  S  (yellow)  far  away 
from  the  two  spheres  that  corresponds  to  the  scattered  light  inten¬ 
sity  measured  in  our  experiments.  £tip  decreases  exponentially  with 
the  distance,  because  of  the  exponential  decay  of  polariton  modes 
at  the  sample  surface.  However,  Eint  and  S  show  maxima  at  certain 
distances  between  tip  and  sample  of  70  nm  and  50  nm,  respectively. 
The  polariton  modes  at  the  interfaces  of  the  superlens  enhance 
the  evanescent  waves  arising  from  the  object.  Hence,  a  larger  field 
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Figure  6  |  Simulations  of  electric  field  distributions  and  transfer  function  for  the  asymmetric  superlens,  (a)  Planar  superlens  with  a  line  source  as 
object  on  one  side  for  different  wavelengths,  showing  a  confined  field  on  the  image  side  of  the  lens  for  X  =  13.5-14  pm  due  to  superlensing  (all  figures  with 
same  colour  scale  in  arbitrary  units),  (b)  The  isothermal  contour  of  transfer  function  in  the  wavelength  range  of  our  interest  plotted  versus  wavelength  A 
and  wavevector  kv  The  colour  represents  the  transfer  function  (the  square  of  the  ratio  between  the  transmitted  electric  field  after  the  superlens  and  the 
incident  field).  The  white  line  is  the  light  line  in  air.  (c)  Transfer  functions  |T|2  for  the  asymmetric  superlens  (blue)  and  the  control  sample  (red)  at  13.5pm 
wavelength.  The  control  sample  replaces  the  400  nm  SrTi03  film  in  the  superlens  by  a  400  nm  BiFeOs  layer.  One  can  clearly  see  that  the  evanescent  wave 
is  enhanced  by  the  superlens  over  a  large  range  of  wave  vectors  (up  to  10  k0\  The  sharp  peaks  around  k0  are  due  to  total  internal  reflection,  (d)  A  planar 
superlens  with  two  spherical  objects  on  both  sides  for  A  =  14  jim  and  increasing  gap  z  between  the  upper  sphere  (probe)  and  the  sample  surface,  (e) 
Parameters  of  interest  extracted  from  simulations  as  shown  in  d:  the  electric  field  at  the  lower  apex  of  the  probe  Etlp,  the  electric  field  at  the  SrTi03-BiFe03 
interface  Eint  and  the  integrated  Poynting  vector  S  far  away  from  the  two-sphere  system.  All  parameters  are  depicted  as  a  function  of  the  gap  z.  In  contrast 
to  Etip,  which  has  the  highest  value  at  z— >0,  Eint  and  S  show  maxima  for  certain  z  being  70  and  50  nm,  respectively. 


amplitude  of  polariton  modes  Eint  for  a  certain  probe-sample  dis¬ 
tance  indicates  enhanced  superlensing  of  the  coupled  probe-object 
system  for  this  position  of  the  probe  scatterer.  Moreover,  the  scat¬ 
tered  light  (represented  by  S)  from  this  system  has  a  maximum 
around  the  same  gap-size,  which  correlates  the  far-held  observation 
with  the  enhanced  coupling. 

Discussion 

The  theory  of  transformation  optics46-48  states  that  a  superlens  with 
negative  index  of  refraction  distorts  the  optical  space  in  a  remark¬ 
able  way:  the  space  is  folded  by  the  superlens,  with  the  object  plane 
and  the  image  plane  (as  well  as  a  plane  within  the  lens)  at  the  same 
position  in  optical  space46.  If  we  place  scatterers  in  both  planes,  one 
can  consider  them  to  be  two  induced  dipoles  located  at  exactly  the 
same  position  in  optical  space,  which  therefore  might  lead  to  the 
observed  enhanced  coupling  of  both  scatterers. 

We  note  that  a  maximum  in  the  image  plane  of  a  superlens  was  also 
observed  earlier  in  the  microwave  range  using  loop  antennas  as  source 
and  detector49-51.  The  position  of  the  maximum  changed  because  of 
matched  coupling  between  object  and  probe,  when  using  detector 
antennas  with  different  radii  and  different  resistance  load50.  However, 
any  dependences  on  the  wavelength  and  on  corresponding  changes  in 


the  sample  dielectric  constants  were  not  studied.  It  was  proposed  that 
this  effect  could  be  applicable  for  three-dimensional  imaging50,  which 
might  be  possible  with  the  oxide-based  superlens  as  well. 

In  this  paper,  we  proposed,  designed  and  demonstrated  per- 
ovskite-based  superlenses  for  electric  fields.  These  materials  are  in 
particular  suitable  for  superlenses  in  the  infrared  range,  showing 
matching  pairs  of  the  real  parts  of  their  dielectric  constants  with 
opposite  signs.  Moreover,  perovskites  show  low  intrinsic  absorption 
at  the  wavelengths  of  interest  and  different  materials  can  be  grown 
epitaxially  on  top  of  each  other.  Both  effects  lead  to  low  losses  in 
the  superlens  structure.  As  some  perovskite  oxides  are  ferroelectric 
(for  example,  BiFe03  in  the  current  study),  a  superlens  consisting  of 
these  materials  might  be  tunable  by  an  external  electric  field. 

We  characterized  the  perovskite  superlenses  in  both  spectral 
and  three-dimensional  spatial  domain.  The  near-field  examination 
shows  sub  wavelength  resolution  of  A/ 14  at  the  superlensing  wave¬ 
length.  Comparison  of  symmetric  and  asymmetric  superlens  show 
a  stronger  superlensed  signal  for  the  asymmetric  case  in  which  the 
probe  is  placed  directly  at  the  superlens  interface.  This  finding  is  sup¬ 
ported  by  the  calculated  transfer  function  of  the  asymmetric  super¬ 
lens,  which  shows  a  better  performance  compared  with  the  symmet¬ 
ric  case  in  terms  of  field  enhancement  and  wave  vector  bandwidth. 
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In  addition,  we  discussed  the  coupling  of  near-field  probe  and 
object,  being  placed  on  two  opposite  sides  of  a  superlens.  It  is  found 
that  the  coupling  effect  is  strongly  enhanced  at  phonon -polariton  res¬ 
onances  and  is  dependent  on  the  probe-object  distance.  These  find¬ 
ings  reflect  the  fact  that  in  a  system  consisting  of  object,  superlens 
and  detecting  probe  coupling  between  all  constituents  takes  place. 
The  observed  superlens-mediated  interaction  between  two  particles 
might  allow  for  controlled  enhanced- coupling  effects.  These  effects 
could  find  potential  applications  in  local  thermal  sensors  as  described 
by  Shen  et  al.52  as  well  as  in  metamaterial-based  multifunctional  cir¬ 
cuits53’54.  We  envision,  for  example,  an  enhanced  transport  of  polari- 
tons  from  particle  to  particle  by  these  superlenses:  a  series  of  super¬ 
lens-coupled  particles  can  act  as  a  polariton  conductor  with  narrow 
bandwidth,  which  is  determined  by  the  superlensing  material. 

Methods 

Growth  and  absorption  losses  of  perovskite  oxides.  The  various  thin  films  were 
grown  by  pulsed  laser  deposition.  Thin  films  of  SrRuOs  were  grown  at  680  °C  and 
films  of  BiFeOs  and  SrTi03  were  grown  at  700  °C  in  lOOmTorr  of  oxygen55.  Follow¬ 
ing  the  growth  of  the  SrRu03  layer,  features  were  defined  using  photo -lithography 
and  samples  were  ion-milled  to  produce  the  objects;  subsequently,  BiFe03  and 
SrTiOs  films  were  grown.  Films  were  found  to  be  single  phase  and  fully  epitaxial  in 
all  cases  via  X-ray  diffraction  and  transmission  electron  microscopy. 

To  compare  the  absorption  of  our  perovskite-based  superlens  with  former 
materials  used,  we  calculate  the  ratios  n/  k  with  n  and  k  being  the  real  and  imagi¬ 
nary  parts  of  the  refractive  index  at  the  wavelengths  for  which  superlensing  is 
expected.  Larger  n/K  corresponds  to  less  material  absorption  losses.  For  SrTiOs 
this  ratio  is  12.7xl0-2,  three  times  larger  than  for  SiC  and  silver,  which  both  show 
a  n/K  of  about  4.2xl0“2  at  their  corresponding  superlensing  wavelengths18,17.  Note 
that  the  superlensing  condition  for  the  perovskite-based  superlenses  is  fulfilled  for 
eSrTio3=_  2,  whereas  for  the  SiC  and  silver  superlens  it  was  observed  around  £=-3. 
However,  even  when  £SrTi03=-  3  at  X-  14.7  pm,  we  find  n/K-  10.5xl0'2  being  much 
larger  than  for  the  other  superlenses. 

s-NSIM  setup.  In  s-NSIM,  a  scattering  probe,  namely  a  metal-coated  atomic- 
force  microscope  tip  with  a  typical  radius  of  50  nm,  is  placed  close  to  the  sample 
surface34,35.  This  probe  transforms  the  evanescent  fields  into  propagating  waves, 
which  can  be  detected  in  the  far-held34,56.  To  separate  the  near-held  from  the 
much  larger  far-held  signal,  we  use  the  method  of  higher-harmonic  demodula¬ 
tion56,57:  in  tapping-mode  atomic-force  microscope,  the  distance  between  tip  and 
sample  is  modulated  with  small  amplitudes  of  about  30  nm  at  frequencies  around 
£2=  150  kHz.  On  this  scale,  the  far- held  changes  linearly  with  distance  resulting  in 
a  modulation  with  the  same  frequency  Q.  On  the  other  hand,  evanescent  waves 
depend  non-linearly  on  the  distance,  resulting  in  a  modulation  at  nQ  ( n  =  1,2,3,. . .). 
Hence,  when  hltering  the  modulated  signal  at  n£l  with  n> 2,  at  higher  harmonics, 
we  obtain  near-held  components  only,  with  less  far-held  contributions  for  larger 
orders  of  n.  If  not  stated  otherwise,  the  near- held  signals  shown  in  this  paper 
represent  third-harmonic  signals  NF3Q  to  ensure  pure  near-held  examination.  For 
comparisons  at  different  wavelengths,  we  nomalize  NF3n  to  the  current  laser  power 
and  to  the  spectral  response  of  detector  and  optical  elements  in  the  beam  path. 

Our  investigation  was  performed  using  the  free-electron  laser  FELBE  at  the 
Helmholtz-Zentrum  Dresden-Rossendorf  (www.hzdr.de),  Germany,  which  offers 
continuous  tunability  across  a  wavelength  range  of  4-250  pm  at  an  average  power 
of  up  to  10  W  (delivered  as  a  picosecond  pulse  train  at  a  repetition  rate  of  13  MHz). 
In  the  10-20  pm  wavelength  regime,  the  typical  spectral  full  width  at  half 
maximum  (FWHM)  of  the  FEL  is  50-100  nm. 

Near-field  interaction  and  probe-sample  coupling.  When  imaging  a  superlensed 
signal  by  s-NSIM,  the  near-held  probe  can  not  be  assumed  to  be  a  passive  element 
only,  but  the  coupling  between  probe  and  objects  needs  to  be  taken  into  account. 
The  origin  of  the  coupling  lies  in  the  optical  interaction  principle  of  NSIM,  which 
results  in  the  consequences  for  the  NSIM  examination  of  superlenses  as  discussed 
in  the  following. 

In  our  s-NSIM  and  the  similar  one  used  in  ref.  18,  the  probe-sample  system  is  illu¬ 
minated  from  the  probe  side  and,  hence,  the  initial  near  held  is  generated  by  the  probe 
rather  than  by  the  sample.  The  function  of  the  probe  is  two-fold:  hrst,  the  evanescent 
waves  arising  from  the  probe  illuminate  the  sample  locally  and  excite  modes  in  the 
sample  and  on  its  surface.  Second,  the  probe  senses  the  helds  of  the  sample  modes 
and  transforms  them  by  scattering  into  propagating  waves,  which  can  be  detected  in 
the  far-held.  This  probe-sample  interaction  depends  on  the  scattering  behaviour  of 
the  probe,  the  local  optical  properties  of  the  sample,  as  well  as  the  distance  between 
the  probe  and  the  sample  surface.  As  a  hrst-order  approximation,  it  can  be  described 
by  the  dipole  model  introduced  by  Knoll  and  Keilmann56:  here,  the  s-NSIM  signal 
is  described  by  scattering  in  the  near-held  coupled  probe-object  system  represented 
through  an  effective  polarizability.  A  more  sophisticated  approach  describes  the  probe- 
sample  interaction  through  higher-order  modes58  showing  the  same  qualitative  results. 


A  critical  point  for  successful  NSIM  measurements  is  that  the  probe  proper¬ 
ties  and  illumination  are  kept  constant  during  the  measurements,  therefore  the 
changes  in  the  near-held  signal  are  essentially  caused  by  local  sample  properties.  In 
our  study,  we  particularly  keep  the  position  of  the  probe  hxed  while  the  sample  is 
scanned  to  ensure  constant  illumination.  The  probe  properties  like  material,  geom¬ 
etry  and  orientation  with  respect  to  the  polarization  of  the  incident  light  are  hxed 
as  well.  Hence,  the  observed  contrast  arises  from  changes  in  the  sample  only.  The 
inhuence  of  probe  or  object  properties  are  enhanced,  when  either  one  of  the  two 
is  excited  close  to  its  resonances.  In  the  mid-infrared,  a  metallic  probe  shows  no 
distinct  resonances,  whereas  the  samples  discussed  in  our  work  support  polariton 
resonances.  Hence,  it  is  the  sample  that  determines  the  NSIM  signal,  whereas  the 
tip  acts  solely  as  a  small  optical  dipole. 

In  the  case  of  a  superlens  as  the  sample,  the  near-held  interaction  is  extended: 
at  the  superlensing  wavelength,  the  evanescent  waves  arising  from  the  probe  are 
transferred  to  the  object  plane.  The  objects  are  excited  by  these  waves  and  generate 
additional  evanescent  helds,  which  are  reconstructed  by  the  superlens  on  the  image 
side  of  the  sample.  Finally,  the  near-held  signal  at  the  image  plane  is  transformed 
by  the  probe  into  detectable  propagating  waves.  As  the  observed  signal  depends 
strongly  on  the  properties  of  the  superlens  as  well  as  the  positions  of  both  scat¬ 
tered,  probe  and  object,  with  respect  to  the  superlens  slab,  we  call  this  effect  a 
superlens-enhanced  probe-object  interaction. 

Besides  the  superlens -enhanced  interaction  between  probe  and  objects,  the 
probe  can  also  interact  with  the  surface  layer  of  the  superlens,  which  is  in  our  case 
either  BiFeOs  or  SrTiOs  for  the  symmetric  or  asymmetric  superlens,  respectively. 
We  observe  enhanced  near-held  signals  due  to  propagating  polariton  modes  at 
the  toplayer-air  interface  of  the  samples  at  wavelengths  around  17  and  14.6  pm 
for  symmetric  and  asymmetric  superlens,  respectively.  Please  note  that,  unlike 
the  localized  modes  at  the  superlensing  wavelength,  this  signal  does  not  carry  any 
information  about  the  SrRuOs  objects  (see  Figs  3  and  4). 

Numerical  simulations.  The  numerical  simulation  results  presented  in  the  manu¬ 
script  are  all  based  on  the  commercial  hnite-element  solver  COMSOL  3.5.  Because 
of  memory  limitation,  the  simulation  is  performed  in  two  dimensions.  The  extrem¬ 
ities  of  the  simulation  domain  are  assigned  scattering  properties,  which  essentially 
mimic  the  necessary  open  boundary  conditions.  We  modify  the  size  of  simulation 
domain  and  (local)  meshes  to  ensure  that  the  held  variation  is  less  than  1%.  Such 
a  convergence  of  the  numerical  simulation  verihes  that  the  boundary  condition 
and  meshing  are  assigned  properly.  The  total  mesh  number  is  up  to  450,000.  For  a 
workstation  with  32G  RAM  memory  and  two  dual  CPUs  (2.66  GHz),  each  simula¬ 
tion  takes  less  than  2  min  to  converge  and  reach  the  relative  tolerance  of  10'6. 

To  numerically  retrieve  the  s-NSIM  signal  shown  in  Figure  5,  iterative  simula¬ 
tions  with  varying  tip  positions  are  performed  by  combining  a  Matlab  script  with 
COMSOL.  The  near-held  probe  is  made  by  gold,  whose  dielectric  constant  is  taken 
by  htting  the  data  from  Johnson  and  Christy59.  The  probe  is  modelled  as  a  triangle, 
whose  tip  angle  is  about  22°  and  tip  apex  is  a  half  sphere  with  the  diameter  of 
100  nm.  The  superlens  consists  of  layers  of  BiFeOs  and  SrTi03  with  thicknesses  of 
200  and  400  nm,  respectively,  and  the  dielectric  constants  are  taken  from  litera¬ 
ture20,21.  The  SrRu03  objects  are  assumed  to  be  rectangles  of  50nmx4pm  with  a  die¬ 
lectric  constant  of  SrRu03  determined  by  Fourier  transform  infrared  spectroscopy 
(see  Supplementary  Fig.  SI).  The  system  is  illuminated  by  monochromatic  light  with 
an  incident  angle  of  75°.  To  study  the  response  of  the  system,  we  calculate  the  elec- 
tric-held  distribution  at  all  positions  of  the  system.  We  observe  a  held  enhancement 
underneath  the  tip  and  at  the  superlens  interfaces  because  of  polariton  excitation. 
The  scattering  cross-sections  (Fig.  5b)  as  measured  in  the  experiment  are  related  to 
the  integrated  Poynting  vector  over  a  solid  angle  about  70°  from  the  tip  apex.  Such 
a  simulation  is  repeated  for  different  tip  positions,  in  direct  analogy  to  the  scanning 
process.  To  simulate  the  higher-harmonic  demodulation,  we  calculate  the  gradient 
of  the  Poynting  vector  to  obtain  the  hrst-harmonic  signal.  Consequently,  the  gradi¬ 
ent  of  the  hrst  harmonic  reflects  the  second-harmonic  signal  and  so  on60. 
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